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^ I This paper reports a detailed description of the equivalent linear 

VO ■ two-body method for the many body problem, which is based on 

an approximate reduction of the many-body Schrodinger equation 
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lO ' method it has been applied to the one-dimensional A^-body prob- 
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\0 • lem with pair-wise contact interactions (McGurie-Yang A^-body 
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0\ \ problem) and to the dilute Bose-Einstein condensation (BEG) of 

atoms in harmonic traps at zero temperature. For both cases, it 
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^ i is shown that the method gives excellent results for large A^. 
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I. Introduction 

In this paper we present an approximate metliod of obtaining the eigen- 
value solutions of the system of interacting A^ bosons using an equivalent 
two-body method similar to that used by Feshbach and Rubinov [1] for the 
triton {^H) three-body (A^=3) bound state. They [1] used both the varia- 
tional principle and a reduced coordinate variable ( not the hyperradius) to 
obtain an equivalent two-body equation for the three-body bound state {^H). 
For many-body problems, use of one reduced coordinate variable (the hyper- 
radius [2]) was made to obtain equivalent two-body equations by keeping 
only a finite sum of terms of the hyperspherical expansion with K = Kmin 
{K is the global angular momentum). This method has been applied to the 
ground state of the A^-body system composed of distinguishable particles 
or of bosons and also to nuclear bound states [3,4]. It was shown that the 
method leads to the correctly behaved nuclear bound states in the limit of 
large A {A is the nucleon number) [4] . Recently, it has been used to describe 
the Bose-Einstein condensation (BEC) of atoms in isotropic harmonic traps 
[5]. We note that Morse and Feshbah [6] have used hyperspherical coordi- 
nates for solving the problem of two interacting particles in a central field of 
force. More references can be found in [7,8]. 

For the A^-body system, our method for obtaining the equivalent linear 
two-body (ELTB) equation consists of two steps. The first step is to give the 
A^-body wave function ^(ri, r2, ...) a particular functional form 

vI/(ri,r2,...)^'^(Ci,C2,C3), (1) 

where (1,(2,^''^^ Cs ^ire known functions. We limit C's to three variables in 
order to obtain the ELTB equation. The second step is to derive an equation 



for ^(Ci, C2) Cs) by requiring that \I/ must satisfy a variational principle 

6f^*^dT = (2) 



with a subsidiary condition /^*^(ir = 1. This leads to a linear two-body 
equation from which both eigenvalues and eigenfunctions can be obtained. 
The lowest eigenvalue is an upper bound of the lowest eigenvalue of the 
original N-body problem. 

In Section II, we apply the method to the one-dimensional A^-body prob- 
lem with pair-wise contact interactions (the McGuire-Yang A^-body problem 
[9, 10]) and demonstrate that our method is a very good approximation for 
the case of large A^. In Section III, we consider the Bose-Einstein condensa- 
tion (BEG) of atoms in isotropic harmonic traps. In Sections IV and V we 
apply the method to the dilute BEG of atoms in anisotropic harmonic traps. 
It was shown the method gives excellent results for large A^. In Section VI 
we investigate the stability of the BEG for the case of atoms with negative 
scattering length at zero temperature using our method. A summary and 
conclusions are given in Sec. VII. 

II. McGuire-Yang N-body Problem 

There are only several known cases of exactly solvable three-body and 
four-body problems. For the A^ = 3 case it was shown [11] that the Faddeev 
equations [12] for one-dimensional three-body problem with pair-wise con- 
tact interactions are exactly solvable. For the one-dimensional A^ = 4 case, 
analytical solutions of the four-body Faddeev-Yakubovsky were obtained in 
[13]. We note that for nuclear three-body systems with short-range inter- 
actions, the Schrodinger equation in three dimensions is reformulated into 



the Faddeev equations [12] which have been solved numerically after making 
partial wave expansion [14] or without partial wave expansion [15]. In the 
following, we consider an exactly solvable one-dimensional A^-body system 
as a test case for our method. 

For the one-dimensional N-body problem with the Hamiltonian 

1 ^ 

i=l i<j 

the Schrodinger equation 

H^ = E^ 

is exactly solvable. The bound and scattering states for this system have 
been found by McGuire [9] and by Yang [10]. 

For the case c < 0, there are bound states [9] for the system of N bosons 
with the wave function of the following form 



vi> = exp[(c/2)^ 






J 



and the energy of this bound state is given by 

E = -c^N{N^ - l)/24. (3) 

The McGuire- Yang (MY) N-body problem provides a unique possibility 
of checking the validity and accuracy of various approximations made for the 
Schrodinger equation describing A^ particles interacting via short range po- 
tential. In the following, we describe the equivalent linear two-body equation 
(ELTBE) method and compare its solution with the exact solution, Eq. (3), 
of the MY problem. 

For this case, we seek for eigenfunction ^ of iJ in the form of 



where 

P 



\ 



1^ 2^\^i ^j) ■ 

i<j 



We now derive an ELTB equation for ^(p) by requiring that \E' must satisfy 
a variational principle (2). This requirement leads to the equation 

h^ + ^^^^ + 2V^(P)]^(P) = 2i5^(p), (4) 

Vip) = ~9/p, (5) 



where 



, r A^/2 - 1/2 

g = cN(N - 1 ^' '—^ 6) 

^ ^2V2^r(A^/2-l) ^ ^ 



with 



(see Appendix A for details). Eq. (4) with the Coulomb like potential V{p) = 
g/p, Eq. (5), can be solved analytically. 

Instead of the parameter E and variable p in the Schrodinger equation 
(4), we introduce the new quantities 



^=JP, V = 9^^ (7) 

For negative g and negative energies, f] is real negative number. On making 
the substitutions (7), Eq. (4) becomes 

To calculate the asymptotic behavior of ^ for large z, we omit from Eq. (8) 
the terms in 1/z, and l/z"^, and obtain the equation 

dz^ 4 ' 



which shows that ^ behaves as exp[~z/2\^ a correct asymptotic form. 
After substitution 

Eq. (8) becomes 

d^V r^rr. -.dy , N — 2. ^ .. 

,_| + (iV-2-.)-|-h + ^), = 0. (9) 

The solution of this equation (finite for z = 0) is the confluent hypergeometric 

function 

A^ -2 
y = F{r]+^^,N-2,z). (10) 

A solution which satisfies the condition at infinity is obtained only for nega- 
tive integer (or zero) values of r^ + ^^, 

N -2 
r/ + ^— = -n,(n = 0,l,...). (11) 

From the definition, Eq.(7), of the parameter t], we find 

EJN) = -2- -. (12) 

"^ ' {N-2 + 2ny ^ ' 

Using Eq.(6), we have the following expression for ground state energy 

F(m- o\ N{N-l)T{N/2-l/2) 

We note that our system has only one bound state, n = 0, consistent with 
the exact solution, Eq.(3). States with n 7^ are not bound states because 

\EnMN) |<|^o(A^-l) |. 

In the case of large A^, using the asymptotic formulas for F function, 

lta^^ = .°-^(l + 0(i)). 



we obtain 



,2 



Eo{N) = N' 



3 



in 



for the leading term of Eq. (13). On the other hand we have for the large N 
case from Eq. (3), 

E--—N'^(l + 0(^)). 

- 24 ^ ^iV2^^ 

Therefore, for the McGuire-Yang N-body problem, we have demonstrated 
that the ELTBE method, Eqs. (4-6), is a very good approximation for 
the case of large A^ (the relative error for binding energy is about 4.5 %). 
Furthermore, our approximation, Eq. (13), agrees remarkably well with the 
exact value, Eq. (3), for any N (the maximum value of relative error for the 
binding energy occurs for A^ = 3 and is about 10 %). 

III. Isotropic Trap 

In this section, we consider A^ identical bosonic atoms confined in a har- 
monic isotropic trap with the following Hamiltonian 

^2 TV . TV 

H = -—J:^^ + ^rnu^Y^rl + J:VUr^-r,). (14) 

i=l i=l i<j 

For the eigenfunction \1' of H, we assume that the solution for ^ has the 
following form 

^(-"'i'--^)-;^(I^' (15) 



where 



N 

p = Y.rl (16) 



We now derive an ELTB equation for \& by requiring that ^ must satisfy the 
variational principle (2). This requirement leads to the equation 



where 



H^ = E^, (17) 



h^ d^ m 2 2 ^^ (3N-l)(3N-3) ,,, , 



with 



^^ v^ r(3iV/2 - 3/2) p3 io ^ '^ 2p2^ 



(19) 
(see Appendix B for details). 

We note that Eq. (17) is exactly the form of the Schrodinger two-body 
equation in which a centrifugal potential energy is given by (A^ — 1)(A^ — 
3)/(4p2) with identification of angular momentum quantum number / = (A^ — 
l)/2. 

In the dilute condensate case 

P > TA, (20) 

where ta is an atom-atom interaction range. Hence we can use for V{p), 
Eq.(19), the following approximation 

N{N-1) r(3A^/2) 1 f^ 2 



"•'">= VS r(3iV/2l3/2) ;;^/o "•"'<'■>''*■ P^' 



which is proportional to the scattering length a^ in the Born approximation 
for binary collisions 






as = —T2 / VMij-)r'dr. (22) 

47ra Jo 



Since the actual atom-atom interaction in condensate is much larger than the 
scattering energy, it is not possible to use perturbation theory to describe 
the scattering. However, since large changes in the wave function only occur 
over very small distances and since the wave function outside the range of 
interaction is only slightly changed by the interaction, Fermi [16] realized 
that it is possible to introduce a pseudopotential, which can be used to 
calculate small changes in the wave function outside the range of interaction 
by perturbation theory. 

Using the above argument, let us introduce a pseudopotential U, so that 
the two-body Schrodinger equation becomes 

A^ + Uij = E'^, (23) 

m 

with 

{Uo ii r < f], 
(24) 
otherwise, 

where r] is some distance chosen so that rj ^\ a \, rj ^ ta, and rj <^ p, with 
the scattering length a and range of condensate p. 

We hope to find U so that the exact atom-atom wave function for r > rj 
will be given by the solution of Eq. (23) with amplitude f{9) given by the 
Born approximation applied to U 

fsie) = --^ I rf\f/(r)e^*-^^/)-. (25) 

Since we know that the exact f{9) is independent of 9 for our very low energy 
case, the range rj of the pseudopotential U{r) must satisfy 

A;r7 < 1, (26) 



where k = JmE/h . In this case Eq. (25) reduces to 

771 

fsie) = -as = -^Uov'- (27) 

In order for the Born approximation to be valid we require the following 

relation 

< 1, 28) 

a 

where a is the exact scattering length with pseudopotential U{r). 
Combining Eqs. (26-28) we find 

770 < ^ < ^, (29) 

where ?7o is fixed from Eq.(28). 

Therefore, the concept of pseudopotential may be used as long as the 
energy of relative motion, E, is small 

E < -^ (30) 

mriQ 

where r]o is fixed from the condition "~"^ < 0.1. For the ^Li case we have 

rio = 0.2 X lO^^cm, (31) 

and for the case of ^^i?6 atoms we have 

r]o = 0.7 X 10~^cm. (32) 

Eqs. (30-32) imply that the temperature of the condensate must be consider- 
ably lower than lO""^-^' for the ^Li condensate and 7.6 x 10~^K for the ^^Rb 
condensate. 
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IV. Anisotropic Trap 

In this section, we consider N identical bosonic atoms confined in a har- 
monic anisotropic trap with the following Hamiltonian 



^2 JV . N 

H = -7r-T.^^ + 7;i: raiyjlA + ^Ivl + ^l^]) + E ^Un - r,), (33) 

For eigenfunction \E' of H, we assume the solution for \1' has the following 
form 

where 

N N N 

i=l i=l i=l 

We now derive an equation for ^(x, y, z) by requiring that \^(a;, y, z) must 
satisfy the variational principle (2). This requirement leads to the equation 

M = E^, (36) 



where 






(37) 



with 



^ '^' ^ 2(27r)3/2 ^r((iV-l)/2)^ a;|/;2 v ,y, ;, \ j 



11 



and 

G{x,y,z) = I dx' dy' dz'VU^{x'f + {y'f + {z'?) 

J-V^x J-V2y J-V2z ^ 

(39) 
(see Appendix C for details). To the best of our knowledge, Eqs. (38) and 
(39) have not been discussed in the literature. 
In the dilute condensate case 

X :^ ta, y ^ rA, -2 > ta (40) 

we can use for V(x,y,z) the following approximation 

N{N-1) T{N/2) 3 1 



V{x,y,z) 



2(27r)3/2 'r((A^-l)/2)' xyz 



(41) 



/\/2x rV^y (•\/2z I 

dx' ^ dy' / dz'VU^{x'f + {yr + {z'Y). 
-V2x J-V2y J-V2z ^ 

Since we have the following relation for the case of large x, y, and z, 

dx' dy' dz'VU^ix'Y + {y'f + {z'f) = 47r / V,nt{r)r^dr 

(42) 

V(x,y,z) is proportional to the scattering length in the Born approximation 

for binary collisions 

y(x,|/,z) = ^-, (43) 

xyz 



with 

aBh^N{N-\), r(Ar/2) 



27rm T((Ar-l)/2)' ' 



(44) 



12 



Therefore, it is reasonable to replace expression (44) with a corresponding 
expression proportional to the exact scattering length a for binary collisions 
(Landau replacement [17]). 

ah-'NiN-l) r(iV/2) 
' = V2^m ^r((iV-l)/2)) • ('') 



This approximation is equivalent to the following approximation for Vint 

-,7 1^ ^\ 47r/i^a ^ , . 

V^nt(n - fj) = d{ri - r-j), (46) 



m 

which is the Fermi pseudopotential [16]. 

For the positive scattering length case, a > 0, we look for the solution of 
Eq. (36) of the form 

^{x,y,z) = ^Q,,#(a;)$f (|/)$f (z), (47) 

where Cijk are solutions of the following equations 

/ y -'^ ij k ,lmn^lmn -'-' / j ■^ijk,lmn^lmn V / 

Lm,n l,m,n 



with 



and 



Using 



and 



H.mmn =< $f)$f $f [ H \ $«$(^)$f >, (49) 

A,.,„„=<$f)$f$f [$r^$l?<^)>. (50) 

$«(a;) = x^^-^y^ exp[-mCu{x/a^f/{2h)], 

$f (^) = z(^-i)/2exp[-mcD(^/7fc)V(2^)], 
13 



and 



we have 

OiiaiPjPrnlkln 

(51) 

with ^ = ^^n, n = 2^jLjm/{2Trh)Na, tb = {uj^^ujyUJz)^ a^ = uj.^/u), ay = 
ujy/uj, and a^ = uJz/Cj- 

For the case of large N, \ijk,imn reduces to the 5-function 

and hence 

Hijk,lmn ~ Ediidjrn^kn- (54) 

Using Eq. (54) we have for the ground state energy 

where parameters a,/9, and 7 are solutions of the following equations 

dE dE dE 

For the case of large N we can neglect the kinetic energy term in Eq. (55) 
J5 = — If"= + f,'3^ + f0'^ + 9^1 (57) 
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Substitution Eq. (57) into Eq. (56) gives 



2 2 2^2^ 2^2 2v^^ 2 2v^^ 



"x" = ^^— , "j;/5 = ^^— , "^7 = ^^— • (58) 



Solutions of these equations 



(2v^^a,aj,a,)i/5 (2v^^a,aj,a,)V5 (2v^^a,a„a,^V5 

, P = , 7 = 

(Xa: Oiii Oiz 



•^x^-^y-^zj 



give for the ground state energy 

E 5.2 
Nhu 4 ^ ^ 

We note that Eq.(59) is the exact ground state solution of Eq.(36) for large 
N. For the case of large A^ we can obtain an essentially exact expression 
for the ground state energy by neglecting the kinetic energy term in the 
Ginzburg-Pitaevskii-Gross (GPG) equation [19] (the Thomas-Fermi approx- 
imation [18]) as 

Etf 5 15 ^.2.2 

iV^=7^y^)^"^ (60) 

Comparing Eq. (59) with Eq. (60), we can see that for the case of large 
A^, the ELTBE method is a very good approximation, with a relative error 
of about 8% for the binding energy. 



V. Large N Limit. 

After we have obtained Eq. (37) the next step is to make a proper choice 
for the large N limit of the Hamiltonian. To do this let us rescale variables 

X, y, z 

X = N^/^x, y = N^/^y, z = N^/^z. (61) 
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We can rewrite Eq. (36) as 

r n'^ ( 92 ^2 92 N I m/, ,2-2 _|_ , ,2-2 _|_ , ,2~2\ 

n^ (jv-i)(JV-3) / 1 , j_ , j_\ , an^jv(jv-i) / r(Jv/2) n3 1 i,|, ^ ^^ 

"•"2778 4Ar2 V£2 "T j^2 -r 52; -r ^;vv^ Vr((Af-l)/2)Afi/2 J £pJ^ Af ^ • 

(62) 

In the large N limit, ^ ~ jv2 ~ i^ '^f ^'^^ order of unity and { y{(n-i)I2)n^/'^ )'^ 
is of the order of (1/2)^/^, and Eq. (62) simplifies to 

I- 2;^'i& + i| + Jf' + '''"<*■ «■ ^>i* = !*• («') 

where 

Equation (63) describes the motion of a particle with an effective mass mN"^ 
in an effective potential Veff{x,y,z). Therefore when A^ -^ 00, the effec- 
tive mass of the particle becomes infinitely large and then the particle may 
be assumed to remain essentially stationary at the absolute minimum of 
Veff{x,y,z). Quantum fluctuations are unimportant in this limit and the 
most significant contribution to the ground state energy is given by 

E = NVeffixo, yo, zo), (65) 

where xq, yo, zq are to be obtained from 

dVeff{xo,yo,Zo) _ dVeff{xo,yo,zo) _ dVeffixo,yo, zq) 



a.. %„ a.„ »■ («' 

Obviously Eq. (65) fails if the effective potential does not possess a minimum. 
Instead of variables x, y, z we introduce the new quantities 



a:, = J— 5, y. = W— y, ., = W— 5. (67) 
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On making the substitutions (67),Eqs. (64) and (66) become 

Veff{xt,yt,zt) = -^[{alxlWyy^t+ayt)+-{-^+-^+^)+^——i (6J 
^ ^ "f-t yt ^t ^ xtytzt 



with 



and 



2 1 

2«.^* + 2x1 ~~ 
HV^ + 2,3 = 

2alzt + 3 = 

Z,Zi 


n 


2 ' 

^t ytzt 

n 


yhtzt ' 

h 


zhtyt ' 



(69) 



In the case of large n = 2Jujm/27rhNa we can neglect j(^ + ^ + ^). In 
this case, solutions of Eq. (69) 

^2/5 .f^2/5 .7^2/5 

2 "' 2 "' 2 "' ('TAN 

2a^ 2a^ 2af 

give for the ground state energy, Eq. (65) 

E/{Nhuj) = jh^'\ 

which is identical to (59). Hence we show the semiclassical nature of the 
large N approximation (59). Corrections to the result of the large A^ limit 
for the finite A^ case may be obtained by incorporating in the theory the 
quantum fluctuations around the classical minimum [20]. 



VI. Stability of BEC 

When the scattering length is negative, the effective interaction between 
atoms is attractive. It has been claimed that the BEC in free space is im- 
possible [21] because the attraction makes the system tend to an ever dense 
phase. For '^ Li, the s-wave scattering length is a = (—14.5 ± 0.4) A [22]. For 

17 



bosons trapped in an external potential there may exist a metastable BEC 
state with a number of atoms below the critical value Ncr [23-31]. 
For the a < case, we can see that potential energy in Eq. (37), 

W{x,y,z) = —{uj^x +UJ y +uj^z )4 / , , ^ 



2 ^ "" 2m 4 x^ y^ z^ xyz 

(71) 

for A^ < Ncr has a single metastable minimum which leads to the metastable 

BEC state. We note that for the case of large Ncr, the ELTBE method 

leads to the same Ncr as the variational GPG stationary theory [30]. To 

show this, let us consider an anisotropic trap, u^ = oUy = u±, u^ = \uj±. 

Local minimum conditions A > 0, where A is a matrix with matrix elements 

Aij = d'^W/dxidxj, can be written for this case as 

nV2(5l(5l - n - \^ 6 ±6l / 32 + O(^) < 0, (72) 

where 6z = {2muj±/hNcr)z'^,6± = {2mu±/hNcr)x^, and 

n = 2{muj^/2Txhf'^Ncr \a\ . (73) 

Setting the left-hand side of Eq. (72) to zero and neglecting 0{j^) terms, we 
obtain the following equations for Ncr 



(74) 



and 

n = 5l/^5l[l + {l + ^^\y\ 

o± 

Eqs. (74) are exactly the same as equations for determining Ncr obtained 
from the variational GPG approach [30]. In this reference [30], it was found 
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that the numerical solution of Eqs. (57) for < A < 1 can be interpolated 



as 

n = e-("+^^') (75) 

with a = 0.490419, (3 = 0.149175. Using Eqs. (74) and (75) we have 

Ncr={^^f' -^r-^ (76) 

uj]_m ^l«l 

For an isotropic trap (A = 1) we obtain from Eqs. (57) 5z = 5± = 5'^^"^ 
and n = 5~^'^0.8 ~ 0.535 which are in agreement with the results of Refs. 
[24, 25]. Taking the experimental values of '^Li trap parameters [32], a;_L/27r 
= 152 Hz, and uOz/I-k = 132 Hz we obtain Ncr= 1456. This value of N^r is 
consistent with theoretical predictions [27-31] and is in agreement with those 
observed in a recent experiment [32]. 

We note that the ELTBE method for a general anisotropic trap can be 
improved using a generalization of the hyperspherical expansion 

M/(ri,...r^)= J2 ^S>.(^'^'^)^£-(^-)^x:(^.)>^£(^.), (77) 



where the hyperspherical harmonics Y^^'^IQ^), y^^(i7y), and Y^'^{Qz) are 

N 92 jv 



eigenfunctions of the angular parts of the Laplace operators J2iLi w^i Z)^ '^ 



and YJi=i '§^) respectively. However, we do not expect a fast convergence of 

i 

the expansion Eq. (77) because of nonuniformity of the convergence of the 
expansion of T.i<j Vintifi — Tj) in x, y, and z. 



VII. Summary and conclusions 

In summary, we have presented a method for obtaining an equivalent 
linear two-body equation from the Schrodinger equation for the system of 
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A^ bosons, using reduced variables and variational principle. To access the 
accuracy of the method it has been applied to the McGuire-Yang N-body 
problem for which the exact solutions are known. Our method gives excellent 
results compared with exact solutions. The method has been applied also to 
the dilute Bose-Einstein condensation in anisotropic harmonic traps at zero 
temperature for both positive and negative scattering length. For large A^, 
our method gives excellent results for all these cases. 
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Appendix A 

In this Appendix, we derive Eqs. (5)-(6). To calculate V{p) for McGuire- 
Yang problem we start from definition 

v{p) = ^^^~ ^^ r dx,... r dx^_, r dtvu^^iv^--'^'^''--'"^'/^. 

2, J —CO J —CO J — CO 

(Al) 
where 

/•oo /-oo /-oo .v^iV-i, 2 2^i 

n= dxi... dxjv^i rfte*E„=i(4-p=^)t. (A2) 

J —oo J—oo J—oo 

Using 

/oo „ 

dxe'""' = {-ity/\^/\ {A.3) 

-oo 

and 

lim/ rft-^-— = r(^) ' (A4) 

^-o7_oo (e-^t)'^ I otherwise, 

we have from Eqs. (A.l) and (A. 2) 

>((A-l)/2) 
and 



^ = ^^P^'-'^^T^T^^^ (^-5) 



For the contact interaction Vint{V^x) = c5{\^x) we have 

^(p) . ,iV(iV - 1) r(y-l)/^) 1 (A7) 

^'^^ ^ ^2v^r(A/2-l)p ^ ^ 
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Appendix B 

In this Appendix, we present an outline for evaluating the effective po- 
tential V{p) for the isotropic case. 
Our starting formula is 

(5.1) 

where 

n = f dfdRdfs...dfN r rfte'('-V2+RV2)ig.(El3-n-^')*. (5.2) 

J J—oo 

Using 

/•oo 1 

/ x''-^exp(~fixndx = - — rfi-^/^Tiu/p), (5.3) 

Jo \p \ 

and Eq. (A. 4) we have from (B.l) and (B.2) 

^^ _ (27rr(3/2))^167rp3^-^ 

r(3A/2) ' ^ ' ' 



and 



^^^ v^ r(3A/2-3/2)p3io "^"^"^^^ 2p2^ ^^"*^'^^- 



(5.5) 

Now we calculate the effective potential V{p) for various potentials Vint- 
Substitution of 

VUr) = A5(r) = ^5(r) (5.6) 



into Eq. (B.5) gives 



..,._ ^(^-l) r(3A/2) A 

^^^ 2(27r)3/2 r(3A/2-3/2)p3- ^ ■'' 



For the case of a square-well potential 



V^.nt(r) = <! - (5.8) 

otherwise, 
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the calculation gives the following result 



, V^o/2 if p < W2, 
j{p) otherwise, 



where 



and pFq is the generalized hypergeometric function. 
For the Coulomb potential Vi„i(r) = a/r we obtain 

^ 2iVr(3iV/2) g^ ^ ^ 

3v^r(3iV/2 - 3/2) p 



For a Gaussian potential Vintij-) = Vqc ^ *" 

Vip) = ^°^^^~^^ $(3A3iV/2;-2/jV), (Rll) 

where $(a, 6, ; x) is the confluent hypergeometric function. For 2/5^p^ ^ 1 
Eq. (B.ll) gives 

^^^^ 2 r(3Ar/2-3/2)(/3p)3' ^ • ^ 



and for 2/3^p^ < 1 



V.,,) . Mi^ ,5,3) 



For the Yukawa potential Vint{r) = Vqc '^^ /r 

^^^ v^r(3iV/2-3/2)p3^6(Ar-l) ' 21 , / > / I , h^ P I ) 



.23JV/4^2-3JV/23(^ _ i)0fp4~3A^/2j(3^/2 - 1, v^pp)r(3iV/2 - 3/2)], 

(5.14) 
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In the A^ = 3 case Eq.(B.14) simplifies 

N{N - l)T{3N/2)Vo 



Vip) 



2ttT{3N/2 - 3/2)/iV 



IM' - eifipf + 30 



:ri5) 



_2g-v^MP(2v/2(/ip)3 + 12(/ip)2 + 15v^/ip + 15)]. 



Appendix C 

In this Appendix, we present an outhne for obtaining the effective anisotropic 
potential. For an anisotropic case we introduce 



X 



N 



\ Z^ ^n; y 
\ n=l 



N 



\ 7 . Vny ^ 
\ n=l 



N 



\ n=l 



(C.l) 



we can then write for effective potential V{x,y,z) the following expression 

A(A-l) f^ , , 



V{x,y,z) 



oo J —oo 



OO /"OO 



dx' / dy' / dz%M{xr + (yr + i^T) 



X 



/OO roo roo roo roo roo /"OO 

dRj, I dRy / dRz / dxs... / dx^ / dy^... / (i|/Ar 
-OO J— OO J —OO J —OO J ~oo •/— OO ^— OO 



X /^ c^.3... /"^ rf.. /"^ rft, /°^ dt, r rft3e^[(^')V^H-«^/^-El3 4-^^]*. 



OO J —OO 



xe 



i[(j/')V2+R^/2+Er=3 S/^-S/']t2 „i[(^')V2+i?^/2+E;^^3 4-^^1*3 



(C.2) 
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where 

/oo roo roo roo roo roo 

dx' / dy' / dz' / dR^ / dRy / dR, 
-OO J —CO J —oo J — oo J — oo J —oo 



roo roo 

X / dxs... I dxN I dy^... I dy^ I dzs... / dz]\j j dti 

-cxo J —OO J —OO J —oo J — oo 



X r dU r rft,e*[("')'/2+«^/2+E!L3-n-'l*ie*[fe')V2+i?^/2+Et3 2'^-2^'l*2 



oo ^— oo 



i[{z'f/'2+Rl/2+Y,^^^,^zl-z^]H 



(as) 

Using (A. 4) and (B.3) we obtain 

n = r(l/2)3^(47r)3a;^~2^^-2^iv-2r(^/2)-3, (^.4) 

and 

2(27r)'^/^ r(A'/2 — 1/2) xyz J-V2x J~V2y J-V2z 

X !/.„,( vWTWTW)[(i - ^)(i - |^)(i - ^)]"/^-^/^ 

(C.5) 
Finally, we quote several examples of V{x, y, z) for various potentials Vint- 
For delta-potential, Vint(f) = \6{r) we obtain 

^^'^'^''^^- 2(27r)3/2 ^r(iV/2-l/2)J xyz' ^^'^^ 

and for the Gaussian potential, Vint{r) = VqC^^ *" we have 

V{x, y, z) = M(^L^$(i/2, Ar/2; -2/3V)$(l/2, iV/2; -2/3V) 



x<l>(l/2,iV/2;-2/3V) 

(C.7) 
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